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Complanart of system of polynomial equations 
Andrey Vlasov [j] 
ITEP, Moscow, Russia 
ABSTRACT 

In this paper we study polynomial maps of vector spaces z% — ► A 1 ^ 12 "" 13 z^z^z^ ■ ■ ■ Zi s and their 
eigenvectors and eigenvalues. The new quantity called complanart is defined. Complanarts determine 
complanarity of solution vectors of systems of polynomial equations. Evaluation of complanart is 
reduced to evaluation of resultants. As in linear case, the pattern of eigenvectors defines the phase 
diagram of associated differential equation z% = A l ^ l2 "' ls z^z^z^ ■ ■ ■ Zi s . Theory of such differential 
Q\ equations arise naturally as extension of Lyapunov's theory of stability for solutions of differential 
equations. The results of this work have a number of potential applications: from solving non-linear 
differential equations and calculating non-linear exponents to taking non-Gaussian integrals. 
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1 Introduction 
1.1 Overview 

This paper is about nondinear algebra [HE], which studies nondinear maps Z{ — > A l ^ l2 '" ls z^z^z^ • • • Z{ s 
and systems of polynomial equations. Nondinear algebra is a direct generalization of linear algebra 
[3]. Developing nondinear algebra can help us to take non-Gaussian integrals, and can convert many 
calculations in physics in exact ones. It can also be applied to theory of stability of differential 
equations. Instead of determinants, the main functions/objects of quantitative nondinear algebra are 
discriminants and resultants [HE]- Resultants determine solvability of system of polynomial equations, 
and discriminant determines the degeneracy of nondinear form. In this paper we consider another, new 
measure of degeneracy of system of polynomial equations, called complanart. Complanart determines, 
whether there are n complanar roots of system, for details see sect |1.2.1[ and sect |2.3[ In sect [3] the 



problem of finding nondinear eigenvectors/eigenvalues is discussed. In sect El the applications of all 
derived methods to the theory of polynomial differential equations are presented. The equations of 
such type arise in some degenerate cases in the theory of stability. In sect [5] all derived techniques and 
methods are illustrated on one example of quadratic map of two variables. 



1.2 Symmetric combinations of the roots 

It is well known (see, for example, [1]), that all polynomial symmetric combinations of the roots of 
system of n — 1 homogeneous equations of n variables in principle can be expressed polynomially in 
the coefficients of these equations. But explicit expressions for the roots themselves, with separate 
expression for each root, exist only in simple cases. "Explicit" means the expression, involving only 
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arithmetic roots and algebraic operations. For example, for one equation of two homogeneous variables 
such expressions exist only for degrees from 1 to 4. So, it is sometimes easier and more convenient to 
study the system without finding all roots, but by studying symmetrical combinations of the roots. 
For example, discriminant of a form determines whether the form is degenerate or not, resultant 
(see sect. 1.5 and [21 H]) of a system of polynomial equations determines whether the system has non- 



trivial solution. One more example of successful application of this approach to higher discriminants of 
polynomials described in [5J. Known difficulties of this way of analysis is that there has been developed 
no clear methods of obtaining expression for arbitrary symmetric combinations of the roots through 
coefficients of the system. In this paper new symmetric combination of the roots, called complanart, 
is considered. The evaluation of this quantity is reduced to evaluation of resultants. Calculation 
of elementary symmetric polynomials of roots, namely the generalization of Vieta formulas, is also 
reduced to evaluation of resultants, see |2.1| 



1.2.1 Complanart 

Complanart is a symmetric combination of the roots of a system of polynomial equations. Let 
• • • , f n -i(x) be n — 1 homogeneous polynomials of arbitrary degrees rj of n variables x\, ■ ■ ■ , x n , 
and let A^ 1 ), • • • , are roots of the system: 

h(x) = o 



fn-l(x) 







This system has in general case N = r\T2 ■ ■ - r n _\ roots with at least one non-zero component up to 
overall rescaling, see, for example [I]. Double roots are counted and repeated two times, roots of third 
order - three times etc. Complanart equals: 



C 



N 

n 



Jl 



J 2 



A 



(in) 



•1> 



,=1 



H<«2<"-<«n 



(1) 



Complanart equals iff there is a set of n complanar roots, i. e. there is a set of roots 

A^ 2 \ ■ ■ ■ ,AW without pair of equal indices i\ ^ 12 7^ 13 7^ • • • 7^ i n satisfying 
£ hh—3n\M\M . . . \M = 0. In particular, if there is at least one multiple root, complanart also 
equals zero. In the case n = 2 complanarity of vectors means their collinearity, so complanart reduces 
to ordinary discriminant of polynomial. If the number of distinct roots is less than n, complanart equals 
1. For example, complanart equals 1 if all equations are linear equations. The method of evaluating 



complanart is described in sect 2.2 More details about complanarts can be found in sect 2.3 



1.3 Eigenvectors and eigenvalues 

Eigenvector and eigenvalue of linear maps are known from linear algebra, and have direct analogues 
in non-linear algebra. Non-zero vector Z{ is called eigenvector of A, if it satisfies (with some A): 

2j 1 Zp 2 Zi :i • • • Zi s — A(z)zj 

A, a polynomial of degree s — 1, is called eigenvalue. In non-linear case one can at first find eigenvectors 
and then find eigenvalues by solving linear equations, see |3.2| Since to find eigenvalues one have to 
solve linear on A equations, the set of eigenvalues is a union of planes in the space of all polynomials of 
degree s — 1. This statement can be reformulated using characteristic polynomial of the map, namely: 

Ch A {\) = R{Ai(z) - X(z)zi} 
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C/iyi(A) = iff A is an eigenvalue of A. Characteristic polynomial possesses decomposition on linear 
in the coefficients of A factors, since the set of eigenvalues is a union of planes in the space of all 
polynomials of degree s — 1. This decomposability was firstly stated in [2], but without a proof. The 



proof will be given in 3.3 



Finding eigenvectors is reduced to solving the system of n homogeneous equations of n + 1 variables. 
Such systems possess complanart. Complanart can be used to determine whether the system has 
complanar eigenvectors. Here "complanar" means complanarity of vectors in extended space, i. e. in 



space with additional homogenizing variable, see 3.6 



The number of eigenvectors of non-degenerate map equals c n \ s = s s _ 1 , if there is no degenerations 
such as coinciding eigenvectors or the case when the map is unit map. The formula for c n \ s was stated 



in |2j, but from considerations for diagonal maps. In sect 3.4 this formula is derived in general case. 



1.4 Applications of our approach 
1.4.1 Nonlinear differential equations 

The eigenvectors of a map A l > %2 "' la entirely determine the phase diagram of the system of differential 
equations: 



z- = A tll2 '" %a 



%i\ 



(2) 



If initial condition is proportional to an eigenvector, the solution is very simple, see sect 4.2 



Application in the theory of stability The equations of type ([2J arise naturally when considering 
the stability of the stationary point of system of differential equations. Consider a system of ordinary 
differential equations: 

m = fi(yi,y2,--- ,y n ) 
2/2 = h{yi,V2,--- ,2/n) 



y'n = fn (2/1,2/2, 



,(°) 



,(o) 



.2/nJ 



be a stationary point of this system: fi(y\ 

(0) 



(0) 



v {0) ) 
, yn ) 



o, • • • , fn(y[ 



(0) 



v {0) ) 



0. 



Let y{ ',- ■ ■ ,y. 

Denote now y[ = y% — y)^' ■ To analyse stability of this point and some quantitative properties of this 
stability/instability one can expand /i, • • • , f n around the stationary point: 



Ui 



dfi. 



+ 



1 a 2 /, 



0,0, ••• ,0 



!/j ' 2\d yj dy : yi 



(3) 



The first term in expansion 



9yj 



0,0,- ,0 



is linear map, the second \\ Q y .Q yh y'jy'k ^ s homogeneous 



quadratic map, and so on. The case when there is a non-degenerate linear term in this expansion 
is well-known, this is a subject of consideration of Lyapunov theory of stability (see [6J) with its 



Lyapunov's indices equal to eigenvalues of the matrix 
linear term in this expansion vanishes, 



dfi 



. But for some differential equations 



9y 3 



o,o,- >o 



0,0, ■■■ ,0 

0. In such cases it is necessary to consider the 



terms in expansion of higher degrees. If only the term of the lowest degree is considered, one arrives 
to the system of type The physically motivated example of system of differential equations with 
vanishing linear term will be given in sect. 4.4 The discussion of stability /instability of points with 



vanishing linear term is given in sect 4.3 For example, if the resultant of the main non-linear term 
does not equal to zero, the point is unstable in its complex vicinity. If main non-linear term has a real 
eigenvector, this point is unstable in its real vicinity. 
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1.4.2 Non-Gaussian integrals 

The formula for Gaussian integral and obtained from it Wick theorem are widely used in modern sci- 
ence. In many cases when it is necessary to calculate something Gaussian integrals are used. Gaussian 
integral is the integral: 

+ ?° / — (a- 1 ) J l Ji 

Z(J) = J e- A ^ +J ^d Xl ■■■dx n = -^e ?~ (4) 

— oo 

For example, Feynman diagram technique uses them. In it this formula is extended from ordinary to 
functional integrating. One is really interested in calculating such quantities: 

/ £ty(^(ziMz 2 ) • • • 0(x fc ))es Km+J^x 

< 0(x!)0(x 2 ) • • • 0(x fc ) >= ^77777^77^ — (5) 

J J=0 

which are called corellators. <j> is a field (or fields), L(4>) is a Lagrangian of this field(s) and J(x)(j>(x) is 
called source of the field. If L{4>) contains only terms, quadratic on (f>, for example L{<p) = di<f>d % (j)— m 2 4> 2 
- a lagrangian for free scalar massive field, these quantities are calculated as follows: 

Z(J) ;e / DMfoMaa) • • • <K*k))e* = _J^L_ -l Z HJ,J) ( 6 ) 

J y/det(L ) 

< (p{Xl)(p{X2) ■■■4>{x k ) >free= ~ 



d a a 



dJ{x\) 9J(x 2 ) dJ(x m ) ^ 'J 
Z\j=0 



J=0 



(7) 



The formula ^ is the direct generalization of Q to the case of functional integrating, and determinant 
in it is so-called functional determinant. Lq 1 is called propagator of the field <j>. This result in another 
form is also called the Wick theorem. Now, to calculate ^ one just expands non-quadratic terms 
in the exponent and calculates only correlators in the free theory (i. e. in the theory with quadratic 
Lagrangian). The quantities ^ in this approach are calculated perturbatively. It is more preferable to 
calculate them non-perturbatively, exactly. To solve this problem, it is necessary to evaluate integrals: 

j e (J^i+A i ^ i x j +B^ Xi x j x k +- ) dxidx2 ... dXn (g) 

in the limit n — > oo (n is the number of Xi). Dots in the exponent substitute parts of the Lagrangian of 
higher degrees. The integrals of type ([§]) are called non-Gaussian integrals. Non-Gaussian integrals are 
also studied in [7]. Now there are no simple methods of evaluating such integrals, and the expressions 
for these integrals were obtained in [S] only for several simple cases. The expressions for these integrals 
depend on the number of variables (unlike the gaussian integral), therefore it is not evident how they 
behave in the limit n — > oo. In [S] it was shown, that discriminants of non-linear forms play an 
important role in the evaluation of the non-Gaussian integrals, e. g. they control singularities of these 
integrals. The possible approach to calculating non-Gaussian integrals is to use some form of canonical 
representation of non-linear form. In general case, non-linear form cannot be brought to diagonal 
representation, but the free parameters of transformations can be used to fix up some coefficients of 
the form. A possible variant of such representation of the maps (AT, M m , . . .) under GL n action is 



presented in sect 3.5 Under the action of SO n the canonical representations of forms Aijk, Ajki m , . . . 
and for maps A^,Ai km , . . . are the same, therefore studying canonical representation of map under 
SO n can help us to evaluate non-Gaussian integrals. 
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1.5 Terms and notations 



1.5.1 Homogeneous and non-homogeneous equations 

Polynomial /(xi, . . . , x n ) in variables xi, . . . ,x n is called homogeneous polynomial, if for any A / : 
f(Xx\, . . . , \x n ) = X d f(xi, . . . , x n ). Non-negative integer d is called the degree of /. Any homogeneous 
polynomial can be made non-homogeneous by dividing it by one of the variables to power d, for 
example, by (x n ) d . After such division, the ratios x\/x n , X2/x n , . . . , x n -i/x n can be taken as new 
variables: y\ = x\/x n ,y2 = xijx n , . . . , y n -\ = Xn-l/ x m so the number of variables was decreased 
by one. The variables yi,...,y n _i are called non-homogeneous variables, and x%,...,x n are called 
homogeneous variables. In this paper we will denote by j-th component of i-th root of a system 
of equations in homogeneous variables, and by A® i-th root of one equation in one non-homogeneous 
variable. For example, if the equation were ao(^2) a! + a\(x2) d ~ 1 xi + • • • + a^-ia^^i) 1 + ad{x\) d = 0, 
then A^ is xi-component of i-th solution, A^ is X2-component of i-th solution. This equation in 



non-homogeneous variable z = f 1 states ao + a\z + • • • + a<i-\z d 1 + adZ d = 0. A w is i-th root 
of this equation. In sect 2 we consider only homogeneous functions of A^ , because all A^ can be 
simultaneously rescaled and remain the solution of the system. Symmetrical combination of roots AW 
is an expression, which does not change under swapping of any pair AW and 



1.5.2 Maps and resultants 

We study homogeneous polynomial maps of vector spaces: 

Zi > *»i Z{ 2 Zi s ■ ■ ■ Zi s 

The degree of the map is denoted by s. If s = 1, we get ordinary linear maps, which are the objects of 
consideration of standard course of linear algebra PJ. Many objects of linear algebra can be generalized 
to describe the non-linear case. First, consider the system of equations 

Aiiiz—is 7 . 7 . „. . . . ~. _ n 

This is a system of n homogeneous equations of n homogeneous variables or n — 1 non-homogeneous 
variables. Such system in general case has no non-trivial solutions at all. Non-trivial solution is a 
solution with at least one component being non-zero. For such solution to exist, the coefficients of 
the system must satisfy one relation, because the number of variables minus the number of equations 
equals one. This relation is 

R{A} = 

where R{A} is polynomial of the coefficients of A, called resultant or hyperdeterminant of A. If s = 1, 
i. e. A is the linear map, the resultant reduces to ordinary determinant of matrix. Resultants play 
an increasing role in modern mathematics and physics. See, for example [21 H El QUI E] for overview, 
[T2"| IT3] for applications in physics and engineering, [14J for application in string theory and [THl HH] 
for computational methods. By non-degenerate map we mean map with non- vanishing resultant. 



2 Complanart and symmetric combinations of the roots 
2.1 Resultant and generalization of Vieta formulas 

Let fi(x),--- , fn-i(x) be n — 1 homogeneous polynomials of arbitrary degrees r« of n variables 
£&!,••• , x n . It is well known, that the system of equations: 

' fi(x) = 

i (9) 
. /n-i(aO = 
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in general case has iV = r\r2---r n projectively-inequivalent solutions, see [I]. Sometimes they may 
coincide, so N should account for multiplicity. For example, one polynomial of 2 homogeneous variables 
or 1 non-homogeneous variab 



e has exactly r\ projectively-inequivalent solutions (with multiplicities). 
Denote different solutions of (9 1 by A* 1 ), A^ 2 ), • 
vector component of i-th solution vector. Thus, Vieta formulas: 



bS N \ Each solution is also a vector, so A^ is a j-th 



symm 



(Ml) A (M2) 



.12 



a (Mjv) 1 
7 Y J ~~ 



A (a(l)) A ( CT (2)) _ _ _ A (<7(iV)) 



.71 



.72 



31— 3 N 



(10) 



Vj 1 ...j N is homogeneous polynomial of coefficients of /«. The degree of Vj v ..j N in the coefficients of 
i-th polynomial fi equals N/r% = r\ ■ ■ ■ rj_irj+i • • • Tj, see [2]. In some particular cases the formula 
for Vj 1 ...j N may be obtained from simple considerations (see for detailed discussion and examples [2]), 
but we consider now a general method of calculating it. This method uses a Poisson product formula 
for resultant, see e. g. [TJ [TB]. Add one more homogeneous polynomial g(x) of degree r of the same 
variables x\, ■ ■ ■ , x n . Then the system of equations 

fi(x) = 



fn-l(x) =0 

g{x) = 

will possess a resultant. Poisson product formula states: 



N 



R{fl,--- ,fn-l,9} =Cl\g(A (i 



ill) 



i=l 



C is a constant, depending on the normalization of the roots; we can normalize the roots in such a way, 
that C will be equal 1. This formula has the following meaning: the system has non-zero solution iff 
g equals zero on one of the roots of other functions. With the help of this formula, one easily obtains 
the tensor V^...^. One can substitute for g linear function: g(x) = g l Xi, then calculate the resultant 



of the system (11) , and (as it is easily seen from Poisson product formula): 



Vji---3N 



the coefficient before g n ■ ■ ■ g jN in R{fi, ■ • • , f n -i,g} 



2.2 Symmetric combinations 



The formula (11) is easily generalized. For example, one needs to calculate: 



N 

P g =U 5(A«, A (J) ) (12) 
At first one calculates the resultant of the system: 

g(x,y) = o 

fi(x) = 
k fn-l(x)=0 

on the variables x\, ■ ■ ■ ,x n , treating y±, ■ ■ ■ ,y n as parameters. This resultant we will denote 

R x {g, fi, ■ ■ ■ , fn-i}- It is still a polynomial in the variables j/i, • • • ,y n - Then one computes the 

resultant of the system: 

Rx{g,h, ■ ■ ■ ,/n-i}(z/) = o 
My) = o 

/n-i(y) = o 
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in variables y\ , 



,y n . Now it is obvious, how to get 

N 



] [ /i(A«,A( J 'U W ) 

i,j,k=l 



N 



i,j,k,m=l 



One should calculate resultant as many times, as there are different arguments in desired function of 



roots. But in sect |2.3| expressions of such a type are needed: 

N 



P' g = []!?(A (! U (J) ) 



(13) 



The difference between (12 1 and (13 1 is, that there is a product over all pairs of non- coincident indices 



in (13), but in (12) there is a product over all pairs of indices. It seems, that (13) can be calculated in 



this way: 



N 



n # (i) ,A w ) 

ns(AW,A(*>) 
1=1 



N 



i,j = l 



(14) 



Numerator in this expression can be evaluated using the formula (12), and the denominator can be 



evaluated using (11). But a problem can arise. If g(x,y) is an antisymmetric function g(x,y) = 
—g(y,x), both numerator and denominator of (14) are equal to zero. So P' g can not be evaluated 



straightforwardly using (14). We have found a way for evaluating PL in this case. Let go(x,y) - some 



antisymmetric linear polynomial on x\, ■ ■ ■ , x n , y±, ■ ■ ■ , y n . Let us consider g(x, y) = go{x, y)+tgi(x, y), 
where gi(x, y) is 



9i(x,y) = ^2 Xi yj = ( ^2 X * J ( 



(15) 



vi=l 



vi=l 



This expression seems strange because it is not invariant under GL n action. We will discuss it a bit 
later. So: 



N 



P;=nSo(A (i) ,AW) = lim 



nS = i«y(AW,Afa')) 

n£i»(AW,A(*)) 



,. nfJ=i(fl0(A (i) ,ACJ-)) + ^(A(O,AW)) 

lim — — 

ni=i(so(AW,A«) + %(AW,AW)) 



(16) 



For non-linear polynomials go(x,y) there is the same technique, the tensor gi(x,y) should be chosen 
in other way. For example, for quadratic go(x,y) : 



9l( x ,y)= ^2 X i X 0VkVm = ( ^2 x i ) ( ^2 x i ) ( ) ( 

i,j,k,m=l 



(17) 



vi=l 



vi=l 



vi=l 



vi=l 
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This method works not only for antisymmetric function of two variables, but for any number of 
variables. This is the formula for three variables: 



N 



N 

P ' 9o ^ II 5o(A«,A (i) ,A W ) = lim. 

i,j,k — l 



N 



II 3(A«,AW,AW) n 0(A«,A«,A«) 

v iJ,fc=l J \i=l 

( Q^(A»,A«,A(i))^ 



Now go(^ 5 - homogeneous function, y, z) = 50(^1 y> 2) + tgi(x, y, z). For linear go: 

n 

gt(x,y,z) = ^2 x iVj z k 

i,j,k=l 



(19) 



These formulas work for g$ of any degree and of any symmetry. If g(x, y, z) is of other degree, g\ 
should be chosen of the same degree. For example, for quadratic go we would write: 



gi(x,y,z)= ^2 XiXjy k y m ziz p 

i,j,k,m,l,p 



(20) 



If the limits (16), (18 1 can be evaluated, they do not depend on the choice of g\. But sometimes, if we 



choose degenerate g\ or simply g\ with small number of non-zero components, both numerator and 



denominator of (16) will be zero even at t ^ 0, and we will not be able to calculate the limit. The 



example of this phenomena will be given in s 2.3.5 Expressions ( 15|17|20 ) are simply examples of 
non- degenerate maps with all components being nonzero. The example of application of the formulae 



(16), (18) is calculation of complanarts. 



2.3 Complanart 

When the system ^ has coincident roots? In the case of two variables (when there is only one 
polynomial) the answer is given by discriminant of the polynomial. Discriminant equals: 



D 



n 



£ km A (i 



ri ( aw - a(j) ) s 



(21) 



i,j=l,i<j 



j,i=l,i<j 



r\ is the degree of a polynomial. The discriminant equals zero iff there is a pair of roots, in which 
one is proportional to another (A^ oc A^,i / j, homogeneous formulation), or two equal non- 
homogeneous roots (AW = \v\ non-homogeneous formulation). In three-dimensional space, however, 
the proportionality of two vectors is defined by two conditions (e. g. X2/X1 = yijy\ and xj,jx\ = yz/vi)- 



But there is a natural one condition for complanarity of three vectors: e 



ijk r 



0. If£ ijk 



%iVjZk 



0, 



three vectors x, y, z lie in one plane, or they are simply complanar. So, we can formulate the condition 
of a system of equations to have three complanar roots: 



N 



C= Yl {£ mlp M 



i,j,k= 1 

i<j<k 



c = o 



(22) 
(23) 



squared for symmetry). The condition itself is (23). (p2J) is some symmetrical polynomial of the 



roots of the system. To make it symmetric combination e mlp Am A^A^ was squared. This symmetric 
polynomial of roots we call complanart. For two variables complanart reduces to ordinary discriminant, 
for four variables: 



C 



N 

n 

i, j,k,q— 1, 

i<j<k<q 



£ mlpr A {i) A U) A (k) A ( q ) 



(24) 



9 



and for n variables: 



C 



N 

n 

ill'" >in=l 
il<«2<---<*n 



e ^2-Jn A (n) A fe)... A ( i „) 
Jl J2 jn 



(25) 



The complanart has degree 2nC^ = 2n ( - jV _^, w , on A, since every factor in product has degree 2n and 



'(N-n) 



n-l 



there are Cjy- factors. Each A oc n (aj) 1//ri , where a« denotes coefficients of /j (see (10)). So: 

i=l 



n-l 



C oc Y\(c 



rsf^n — 1 iV 



(26) 



i. e. deg ai C = 2C]^_ 1 1 ^. If n = 2, then complanart is discriminant, the number of solutions N = r , 
and deg a C = 2 J^J^)! = 2(r — 1). It is well-known expression for degree of discriminant. 

2.3.1 Evaluation of complanarts 

Complanart is just a product of values of antisymmetric function over sets of all different n roots. The 



method of evaluation of such quantities is given in |2.2| It is necessary just to take for go absolutely 
antisymmetric e-tensor of appropriate dimension. The formula for two variables: 



r .. .. nS=i«A»,AW) + %(AW,A(i))) 

C = hm — ^ - - — = Inn 



n&- = ig(A0,Att>) 

t^o n£i^(AW,AW) t^5 n,=i(e(A«,AW) + %(AW,A»)) 



For three variables: 



(27) 



C 6 = lim ■ 



nf, fc=1 p(A« , AO') , AW ) nf=i 5(A» , AW , AW) 



ni-=i5(AW,AW,A0-J 

y> z ) = £ (^ ; y, z ) + tgi(x, y, z) 



(2? 



For four variables: 



N I N 

]1 ^(AW,AW,A( fc ),AM) [ fl 5 (AW,aW,aW,aW) 

,12 ,. i,j,k,m=l \«,J=1 



A? 



C 1 ^ = lim 



n ^(aW.aw.aw.a^; 



JV 



N 



(29) 



I! 5(A«,A»,A0'),A( fc )) Ml ff(AW,AW,AW,AW) 

\i,j,k=l J \i=l / 

g(x, y, z, u) = e(x, y, z, u) + tg 1 {x, y, z, u) 

There is one more fact to explain. Why for three variables the limit yields to C 3 , and for four variables 
it yields for C 12 ? It is because the formulae (27), (28), (29) give us the following expressions: 



A' 



N 



Y[e(A (i) ,A (j) ) instead of ]J (e(A«, A^)) 2 



i+3 



i,3=l 

i<j 



N 



]"| e(AW ) AO') ) AW) instead of JJ (e(A«, A^, A^)) 1 



(30) 
(31) 



i,j,k=l 



i<j<k 



10 



In the case n = 2 these two things coincide (because there is 2 permutations and accounting for them 
adds necessary squaring). In case n = 3 there are n\ = 6 permutations, but a square is again needed - 
so appears C . In case of arbitrary n this procedure yields Ca . 

It is now a simple exercise to write the formula for complanart analogous to ([27]) , (|28|) , ([29]) in any partic- 

N 

ular dimension: it is necessary only that all the factors, appeared in 1 J g(A^ n > , A^ 2 \ • • • , A**™)) 
with some two or more i coincident will be cancelled by factors with explicitly equal i, for example 

N 

etc. 

2.3.2 When complanart is equal to 1? 

Complanart measures linear dependence of n distinct roots of polynomial system of equations fl9J). 
Such system of equations in general case have N = r\r 2 ■ ■ ■ r n _\ projectively-nonequivalent roots, 
where rj., r 2 , • • • , r n _i are degrees of equations. But what happens when A?" < n? It means, that there 
are no n distinct roots, so no n distinct roots can be complanar. Therefore the complanart of such a 
system is equal to non-zero constant, i. e. it does not depend on the coefficients of the system. What 
is this constant? Let us consider the simplest example: one linear equation, for example ax\ + bx 2 = 0. 
It has one solution: A^ = — b, A^ = a. Complanart equals: 



r n£=i<?(A» AW) nfi^AW.AWHtofAW.AW)) 

lim — — = lini — — 

^° n£i<7(A«,A») *-o niI 1 (£(A«,A«) + %(A« 5 A«)) 



But there is only one solution, so the numerator equals g(A^\ A^ 1 )) and denominator equals s^A^, A^ 1 )). 
So complanart equals 1. Due to similar reasons, complanart equals 1 always when 
N — rir2 • • -r n _i < n. For example, complanart equals 1 for any system of appropriate number of 
linear equations. 

2.3.3 Open questions 



The first question is about the choices of g\: (15 1,(171,(20). These formulas are not GL n -invariant. 



The second question: since we do not know appropriate canonical g\, we may take: g(x,y,---) = 
e(x, y,'") + t\gi(x, y, ■ ■ ■ ) + ^252(^5 y, • • •) + ""■ > an d now look at the same limit, but taken on different 



t variables in different consequence. The third question is: by taking a limit at (27 29), only term 
with the lowest degree of t is considered. It is possible, that the terms of higher degrees contain 
information about higher degenerations of the system. For example, vanishing of some higher term(s) 
may correspond to existence of two coinciding roots and so on. This point can be even more interesting 
considering the higher terms in the case of many t. The theory of higher complanarts would also be 
very interesting and useful in applications. 

2.3.4 Examples of complanarts: n = 2, complanarts are discriminants 
Quadratic equation Our first example is: 

f(x) = a(xi) 2 + 6x1X2 + c(x 2 ) 2 , 
or in non-homogeneous variable z = Xi/x2- 

az 2 + bz + c = 

Take g(x, y) = e(x, y) + tg\(x, y) = x\y 2 - x 2 yi + t(xiyi + x\y 2 + x 2 yi + x 2 y2)- For this equation and 
this g\ : 

Y[g(A®,A^) = t 2 {a + c - 6) 2 (4ac - b 2 + t 2 {b 2 - 2ab + a 2 - 26c + 2ac + c 2 )) 

JJ^(AW,AW) = t 2 (a + c-b) 2 
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It is easily seen, that C = 4ac — b 2 , and this expression coincides with discriminant of quadratic 
polynomial. Now take another gi = x±y 2 . For this g%: 

11 g(A®,A®) = t 2 ac(4ac - b 2 + t(4ac - b 2 ) + t 2 (ac - b 2 )) 

hi 

n 5 (A«,A«) = t 2 ac 

i 

The expressions have changed. But the value of complanart is still 4ac — b 2 . 
Cubic polynomial 

f(x) = ax 3 + bx 2 y + cxy 2 + dy 3 (33) 

Complanart of this polynomial is equal 4ac 3 + 4db 3 + 27a 2 d 2 — b 2 c 2 — 18abcd, and it again coincides 
with the discriminant of the polynomial. Other expressions are very long in this case, and we do not 
write them here. 

2.3.5 Examples of complanarts: n = 3 
Two quadratic polynomials, 1 

fi(x) = an(xi) 2 + a 2 ixix2 + ai 3 xix 3 

f 2 {x) = b 22 (x 2 ) 2 + b 2 lX 2 X 1 + b 23 X 2 X 3 

Take, for example, g(x,y,z) = e(x, y, z) + tgi(x, y, z) = {l + t)xiy 2 z 3 -xiy 3 z 2 + {l + t)x 2 y 3 zi-x 2 y 1 z 3 + 



(1 + t)x 3 y\z 2 — x 3 y 2 z\. Than all three terms in (28), namely 

N N N 

I] 3(AW,AW,A W ), 11 sCAW.AWjAW), II f (A W , A W, A W), vanish, so the limit cannot be 

i,j,k=l i=l i,j=l 

evaluated. Nevertheless, if we take, for example: 

n 

9i= ^ Xiy i = XlVl + XlV2 + XlV3 + X2Vl + X2V2 + X2V3 + X3Vl + X3y2 + X3V3 > 



i. e. choose g\ according to (19), this problem is eliminated. The limit (28) equals 



n^=ig(A (t U^,AW)) (n^AW AW, AM) 

t->0 (t-fN 



nS=i5(AW,A«, 

(ai3) 12 (^23) 12 («2i&23 - ai3&22) 6 (fe23«ii - ai36 2 i) 6 (-&22anai3&23 - &2ia2iai3&23 + b 2 ia\ 3 b 22 + bl 3 a 2 ian) 

Complanart can be calculated also by bare hands, i. e. by solving equations fi(x) = 0,f 2 (x) = and 
then substituting these solutions in p2l. Calculated complanart equals: 



(ai 3 ) 4 (^23) 4 (a2i^23 - ai3^22) 2 (ft23an - ai3^2i) 2 (-^22anai3^23 - b 21 a 21 a 13 b 23 + b 21 af 3 b 22 + b 2 23 a 2 iau) 



Thus, the formula (28) holds. 



Two quadratic polynomials, 2 This example is just the particular case of previous one, but it is 
of great importance for further considerations of sect [5] 

fi(x) = x 2 + 2axy — xz (34) 
h(x) = y 2 + 2bxy - yz (35) 



The complanart cubed, and the limit in the formula (J28J) yield: 

C 3 = (1 -2a) 12 (l -2b) 12 (36) 
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3 Eigenvectors, eigenvalues and GL n canonical representation of non- 
linear maps 

3.1 Non-linear eigenvectors, eigenvalues and how they can be found 

Non-zero vector z% is called eigenvector of A, if it satisfies (with some A): 

Af^z h z i2 z i3 ---z is =\{z)zi (37) 

A polynomial X(z) of degree s — 1 is called eigenvalue. Eigenvalue is not a number, because the 
homogeneity constraint is imposed. The homogeneity constraint is nothing but requirement that any 
eigenvector can be multiplied by some non-zero number and remain eigenvector. There is also a 
non-linear analogue of characteristic polynomial of a map, namely: 

Ch A = R{Ai{z) - \{z) Zi } (38) 

By definition of resultant, CHa turns to zero iff exists an eigenvector, corresponding to the polynomial 
X(z). The optimal way to find eigenvectors/eigenvalues in non-linear algebra differs from the way 
in linear algebra. Solving equation CHa = with respect to A is complicated, because it is a non- 
homogeneous equation of M n \ s _i = ^-t)\(s-i)\ coefficients of A of the degree c n \ s = ^5f- It turns 
out, that CHa is not an arbitrary polynomial, it has a structure - namely, this polynomial is always 
decomposable on linear in the coefficients of A factors. It is much simpler to find these factors at 
first, and then find eigenvalues by solving linear equations. To decompose CHai if is necessary to 



find eigenvectors, this is described in sect |3.3| The method of eliminating A from the equation ( |37[ ) is 
described in sect I 



3.2 Zero, nonzero and unitary eigenvectors 

One encounters two cases considering a particular eigenvector zf, X(z) ^ and X(z) = 0. Consider 
the first case. All eigenvectors with X(z) / we call non-zero eigenvectors. This notation can not 
cause any misinterpretation, because each eigenvector is by definition non-zero in the sense that not 
all its components are zeros. Non-zero eigenvector can be rescaled by any non-zero number, because 
the system of equation (37) is homogeneous system. For example, z% can be rescaled as follows: 

Vi = a( Z )vW> - Now > x (y) = A(z Li = !> and 



( A ( 2 )) ' 



My) = Vi 



(39) 



The eigenvectors, obeying (39), are called unitary eigenvectors. A was eliminated from equation for 
eigenvectors/eigenvalues. Thus, all non-zero eigenvectors can be rescaled to unitary vectors, which 

0. Such eigenvectors are called zero-eigenvectors. Any zero eigenvector 



satisfy (39). 

Now consider the case X(z) 
satisfies the equation: 



Ai(z) = 0, 



(40) 



which also does not contain A. 



Any solution of (39) or (40) is an eigenvector. To prove this, it is sufficient to construct A satisfying 
the equation (37) with this z.- L . If yi is the solution of (39), yi is the eigenvector with any A such that: 

KVi) = i. (41) 



as it is easily seen from (37). If z% is the solution of (40), it is eigenvector with any A which satisfies: 

X( Zi ) = (42) 



Thus, the following statement was proven: 
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Proposition 1. All eigenvectors of A can be found among the solutions of (39) or (40), and, otherwise, 



any non-zero solution of (39) or (40) is eigenvector of A 



So, the easy method of finding eigenvector /eigenvalues is: solve equations ( 39p0 ) to find eigenvec- 
tors, and then find A by algorithm, described above. 

3.3 Decomposability of characteristic polynomial 

Since to find A one has to solve linear on A equations (see previous subsection), the set of A, which are 
eigenvalues of A, is a union of planes in the space of all A. Other way to reformulate this statement 
is to say: the characteristic polynomial C7ja(A) is decomposable on linear on A factors. The condition 
of turning to zero one of the factors of CKa defines the plane in the space of all A. These factors are 
the following: 



n-2 



R{Mz) - X(z) Zi } = M(A) H(l- A(e<">)) ]J(X(z^)), 



(43) 



where is the fx-th unitary eigenvector, and n\ is the number of them, z^ u ) is u-th zero eigenvector, 
and ri2 is the number of them, and M(A) depends only on coefficients of A, not on the coefficients of 



A. Indeed, all possible eigenvectors are found among the solutions (39 40) modulo projective rescaling. 



Therefore all possible eigenvalues satisfy either the equation (41) with one of the unitary eigenvectors 



either the equation (42) with one of the zero eigenvectors. If A satisfi es (|4lj ), r. h. s. of (43) turns to 
zero due to one of the factors (1 — X(e fl )). If A satisfies (42), r. h. s. (|43|) turns to zero due to one of 
the factors X(z v ). 

3.4 Number of eigenvectors 



The systems of equations (|39j) and (|40j) can be merged to one system of homogeneous equations by 

(y) s ~ l 



adding an auxiliary variable y: 



Ai(x) 



(44) 



(44) is a system of n homogeneous equations of n + 1 variables x±,X2, • • • ,x n ,y. If a solution of (44) 



has y / 0, we can consider the vector x\/y, ■ ■ ■ ,x n /y, 1, which is a solution both of (44) and of (39), 
and then first n components of it will be an unitary eigenvector. If a solution of ([44 ) has y = 0, the 
first n components x±, . . . , x n of this solution are simultaneously a solution of (40 ) and vector of them 



is a zero eigenvector. The examples of finding eigenvectors using (44) are presented in sect 3.5 and in 
sect [5j So, the following statement holds: 



Proposition 2. All eigenvectors of A can be found among the solutions of (44y (these solutions should 
be normalized, if necessary). 



Consider non-degenerate maps, R{Ai(z)} / 0. Firstly, for such maps the equation (40) has no 
solutions with at least one component being non-zero. The condition R{Ai(z)} / is by definition of 
resultant the condition of non-existence of non-trivial solutions of p0| ). That is, non-degenerate map 
does not have zero eigenvectors. The same statement holds in the linear algebra. This means, that 



n-2 

ri2 = and in (|43|) there is no factor Y[ ^( z v) 

=i 



R{Mz) - X(z) Zi } = M(A) JJ(1 - A(e M )) 



(45) 



Setting in (45) A = 0, one obtains M(A) = R{A}. So 



»i 



R{Ai(z) - X(z) Zi } = R{A} [J(l - A(e (U )) 



(46) 
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For non-degenerate maps, n\ can be easily determined. Consider again (44). If R{A} / 0, it has no 
non-zero solutions with y = (the existence of such solutions would imply the existence of non-zero 



solutions of (40), but this is prohibited by R{A} / 0). As a system of n homogeneous equations (each 
of degree s, s - a degree of ^4) of n+ 1 variables, in general case (44 ) has s n projectively non-equivalent 
solutions, see pi]. It has one solution which is not eigenvector: 



Xi = X2 



y 



(47) 



The number of solutions of (|44|) with at least one x-component being non-zero is s n — 1. Some of these so- 

,(o) 



,x n ). Let 



(0) 



lutions yield projectively equivalent eigenvectors (in the space x\, ■ 

is a solution of (44). Then x±*\ x^\ • • • , Xn\ Us-iU^K where oj s -\ is a root of degree s — 1 from 1, is 
also solution of (44 ). Since in complex plane there is s — 1 distinct roots of 1, s — 1 projectively inequiv- 

x n ,y) solutions of (44) lead to one eigenvector (or in projective equivalent 

This formula was 



alent (in the space x\, 

vectors in x\, ■ ■ ■ , x n ). So the number of eigenvectors in general case is c n i 



previously obtained in [2J, but was obtained from considerations for diagonal maps. So, characteristic 



polynomial has degree c n \ s on the coefficients of A. Some of s n solutions of (44) can coincide (this 
phenomenon is controlled by complanart and higher complanarts, see 2.3). In this case there will be 



less eigenvectors than c n \ s . It should be emphasized one more time, that all counting, leaded to the 
was carried out under the condition R{A} ^ 0. For example, for unit maps (s. 



formula c, 



-1 ' 



3.7.1) all the vectors in the space are eigenvectors. 



3.5 GL n canonical representation of non-linear map 

In linear algebra linear map (and quadratic form) have canonical representation (i. e., become diagonal) 
in the basis of its eigenvectors. The situation in non-linear algebra is similar, but there are some 
differences. Firstly, in general case non-linear map cannot be brought to diagonal from, because the 
group GL n has n 2 parameters, and non- linear map has more parameters. The second difference is, that 



non-linear map in general case has c n \ s eigenvectors, see 3.4 The number of eigenvectors is greater 
than n, the dimension of the space. So there is uncertainty in canonical representation of non-linear 
map, consisting in the freedom to choose one set of n eigenvectors from c n \ s eigenvectors to be basis 
of the space. This uncertainty is yet poorly studied. To obtain canonical representation of non-linear 
map, one should choose as basis vectors any n linear independent zero or unitary eigenvectors of the 



map. The eigenvectors of the system can be found by solving the system of equations (44). Non-zero 
eigenvectors should be normalized to make them unitary eigenvectors. Consider first what happens with 
the components of the map, corresponding to unitary eigenvector. If unitary eigenvector Zj is chosen 
to be a j-th vector of basis, z; t = b\. By definition of unitary eigenvector, Al 112 "" 18 z^z^z^ ■ ■ ■ Z{ s = Z{. 
Substituting z% = 5f, we obtain A\ 3 3 = 5{. If zero eigenvector yi is chosen to be fc-th vector of 
basis, then j/j = 5f. The equation for zero eigenvectors is: A^ n '" li 'y^y^y^ ■■■y% a =0. and therefore 
A\ k '" k = 0. Summarizing the statements about the canonical form: 

AV' 3 = <5| If j corresponds to unitary eigenvector in basis (48) 
A kk "' k = If k corresponds to zero eigenvector in basis (49) 

By this method, only the components with all upper indices equal are fixed, total n 2 components. 
More components can not be fixed, since there are only n 2 free parameters in GL n group. The rest 
components of A remain arbitrary numbers. It is these numbers that determine the map in the non- 
linear case. In linear case the map is defined by n eigenvectors and n eigenvalues. In non-linear case 
the map is defined by n of its eigenvectors (not all!) and by these numbers, namely by the components 
with not all upper indices equal. One more difference in these representations is that in non-linear case 
the normalization of non-zero eigenvectors is important, but in linear case it is not. 
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3.5.1 Example of canonical representation 



The example is the bringing the quadratic map of two variables to canonical representation: 



xi 

■J'2 



A\ l ( Xl ) 2 + 2A\ 2 x lX2 + Af{x 2 ) 2 
A\ l {x l ) 2 + 2A 12 x l x 2 + Af(x 2 ) 2 



(50) 



To obtain canonical representation of non-linear map, one should choose as basis vectors any n linear 
independent zero or unitary eigenvectors of the map. 
obtain the eigenvectors of the map 



As it was established in sect |3.4| and [3^2 
one should solve such systems of equations: 



to 



Af{x x ) 2 + 2Afx l x 2 + Af{x 2 ) 2 
(xx) 2 + 2A l2 xix 2 + Af{x 2 ) 2 



A 11 



yx 2 



(51) 



This is a system of equations (44) for (50), additional homogenizing variable is also called y. Any 
solution of the system (51), (x\,x 2 ,y) with y ^ is non-zero eigenvector for the map (50). Such 
solution can be made unitary eigenvector by dividing all the components of it by y. Then the solution 
will be (x\/y, x 2 /y, 1) and xi/y, x 2 /y are two components of unitary eigenvector. If there is a solution 
of p)I| ) with third component equal 0, e. g. (xi,x 2 ,0), then (xi,x 2 ) are the components of zero 
eigenvector. In some cases, the map can be represented in several forms of (I521 , (|53|) , ([54l . For example, 



if the map has two unitary eigenvectors and one zero eigenvector, it can be represented either in form 



(52) or in the form (53). 



Two unitary eigenvectors are chosen to be basis This choice fixes the components to the 
following values: ^l} 1 = 1, A 22 = 0, A^ 1 = 0, A 2 ? = 1. Canonical representation of the map in 
this case: 



Xl \ _> { (xx) 2 + 2Afx x x 2 
x 2 J ~* V ( x 2) 2 + 2Al 2 Xl x 2 



> i ■) |12 ... ... I ( 52 ) 



One unitary and one zero eigenvector are chosen to be basis vectors In this case "canonical 
basis" consists of one zero eigenvector and one unitary eigenvector. Let unitary eigenvector be the first 
vector of the basis, and zero eigenvector the second. This choice fixes A\^ = 1, A 22 = 0, A^ 1 = 
0, A\? = 0. Canonical representation of the map is: 



Xl \ _> { (xx) 2 + 2Afx x x 2 
x 2 ) ~* \ 2Afx lX2 



12. „ I (53) 



Two zero eigenvectors are chosen to be basis vectors In this case "canonical basis" consists 
of two zero eigenvectors. This choice fixes A^ = 0, A 22 = 0, A^ 1 = 0, A 2 ? = 0. Canonical 
representation of the map is: 



Xl \ _> ( 2A l 2 xix 2 
x 2 J ~* \ 2Afxix 2 



I (54) 



3.6 Application of complanart to the theory of eigenvectors 

As it was explained in |3.3| all the eigenvectors of A can be found by solving this system of equations: 

Mx) = {y) s ~ l Xl (55) 



This is a system of n homogeneous equations of n+ 1 variables, so it possesses complanart (see 2.3). 
Vanishing of complanart of this system means, that there are n + 1 complanar vectors in the space 
(xi,x 2 , . . . , x n , y), i. e. there exists set of n + 1 indices ji, . . . , j n +i with no pair of equal indices, such 
as: 

h...i n+1A (ji) A (h) A 0n+i) = (56 x 
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Here A^^ stands for j-th solution of (55), Aj is i-th component of A^. Here i runs from 1 to n + 1, 

(7) — 
A^ijis y-component of j-th solution. Let one of j indices, for example, ji, corresponds to the solution 

of (55) of type (47), namely {x\ = 0, Xi = 0, . . . , x n = 0, y = 1). Then the sum over i\ in (56) is 

reduced to one term: 



.(n+l)i 2 ...i„ +1 \ji a(J2) a 
IV n+l IV i 2 



,0'n+l) _ Jn+l)i 



, O'n + l) 
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'2,-i'n + l = l 



O'n + l) 



(57) 



The equality A^ + i = 1> following from (47), was used. The explicit symbol of the sum in last line 



is to emphasize, that e now is n-dimensional, not n + 1-dimensional and the sum over indices goes 



,A 



to n, not to n + 1, as in the first line. Expression (57) is up to sign a condition of A 4 - 2 
to be complanar in n-dimensional space of initial coordinates x%, X2, ■ ■ ■ , x n . If the solution (|47]) does 



, (J2) 



(jn+l) 



not enter in the formula (56), the formula (56) does not possess such a simple interpretation. So, 



vanishing of complanart of the system (55) means either that there are n eigenvectors, 
complanar in usual space x\, ... ,x n , or that there are n + 1 eigenvectors, complanar in the 
space with additional homogenizing variable, in space (xi, . . . ,x n ,y). 



3.6.1 Example 



This example is consideration of non-degenerate map of two variables, n = 2 and R{A} / 0. Since 
R{A} ^ 0,the map has only non-zero eigenvectors, which can be rescaled to be unitary eigenvectors. 
Let three unitary eigenvectors (xi,X2,l), (iti, U2, 1), (#i, £2, 1) enter to the formula (56). Then 



formula (56) becomes: 



X1U2Z3 - X1U3Z2 + X2U3Z1 - X2U\Z 3 + - X3U2ZI = X\U2 - X 2 Ul + X 2 Zl - X\Z2 + U\Z 2 - U 2 Zl = 0(58) 



The solution of type (47) equals in this case: (0,0,1). Let (xi, X2, 1), (zi, Z2, 1) and the solution of type 



(47) enter to the formula (56). Then the formula (56) becomes: 



x\Z2 - x 2 zi = 0, (59) 
or simply condition of complanarity of two two-dimensional vectors (x%,X2) and (2/1,2/2) m two- 



dimensional space. Vanishing of complanart of system (55) means, that there are such eigenvectors 
(xi,X2,l), (iti,«2, 1), {z±,Z2, 1), that either (58) or (59) holds. For degenerate map only the equal- 



ity (p8| would change. One more example of using complanart to determine degeneracy of eigenvector 



pattern is given in sect 5.1 



3.7 How to exclude A by one more method 



In sect |3.2| the method of excluding A from equations for eigenvectors by renormalizing vectors was 
explained. There is one more method of excluding A. This method was mentioned in [2]. The method 
is based on following statement: 

Proposition 3. All the eigenvectors of the map A can be found among the solutions of system of 
equations: 



Ai(z)zj - ZiAj(z) = 0, 



(60) 



and otherwise, any solution of (60) is an eigenvector. 



Proof If Zi is an eigenvector of Ai, A{(z) = X(z)zi. Substituting this equality to (60) yields: 
Ai(z)zj — Aj(z)zi = X(z)ziZj — X(z)zjZi = 0. Otherwise, let Z{ is the solution of (60). Multiply both 
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sides of equation (60) by a vector v l , orthogonal to z\. Then one obtains ZjAi{z)v l = Aj{z)ziV l = 0. 
So, Ai(z) is orthogonal to any vector, orthogonal to z%, therefore Ai(z) = C ■ Z{. One can obtain X(z), 
corresponding to this eigenvector by solving the equation X(z) = C. 

At first glance, in (60) there are w ^" 2 ^ non-trivial equations. But only n—1 of them are independent: 



you can fix i and get n — 1 equations for j / i, and all other equations will follow from these. This fact 
was previuosly noted in [2], but without a proof. This is a proposition [I] of sect. 4A and it is proven 
there. 



This is one more method to show linear decomposability of characteristic polynomial (38); one again 



has to solve linear on A equations. For example, for (50) the system of equations (60) reduces to single 
equation: 



A?{ Xl f + {2A? 



Af)x 2 { Xl ) 



+ {Af-2Af){x 2 f Xl -Af{x 2 







(61) 



3.7.1 Unit maps 



Using the system of equations (60), one can establish the condition of each vector of the space to 
be an eigenvector of the map. This condition is condition of vanishing of all the coefficients in the 



equations (60). The map, for which each vector of the space is an eigenvector, is called unit map. Any 
unit map can be represented as Ai(z) = fj,(z)zi, where (j,(z) is a homogeneous polynomial of degree 
s — 1, Characteristic polynomial for unit maps is identically zero, R{A l {z) — \{z)z 1 } = 0. Substituting 
Ai(z) = fJ,(z)zi, one obtains: R{A l (z) — \(z)z 1 } = R{/i(z)zi — \{z)z % ] = R{(fj,(z) — \{z))zi} = 
R{\'(z)zi}. The system of equations \'(z)zi = always has non-trivial solution since \'{z) = always 
has non-trivial solution, because it is one algebraic equation on n variables. So, R{X'(z)zi} = and 
R{A l (z) — X(z)z 1 } = 0. Unit map maps each vector of the space proportional by itself. This is a 
generalization of the notion of unit map in non-linear algebra. For example, the condition of the map 



(51 ) to be unit is: 



2Af 

422 _ 
A 2 



A? 

-< 
2A l * 

Af 



(62) 



These conditions are just the conditions of vanishing of coefficients in (61). 



4 Polynomial differential equations 
4,1 Eigenvectors as stationary points 

As it was mentioned in introduction, when one considers the stability of the stationary point of a 
system of differential equations, one encounters such systems of equations: 

X. /i 1 ^2 y . / . . ... rp . ( f-\ Q ^ 

It is well known, that if the degree of the map equals 1, the solution of this equations is expressed 
through eigenvectors and eigenvalues of linear map A\. In the non-linear case, eigenvectors also play 
an important role in solving these equations, for example, they entirely determine the phase diagram 
of the system. We begin by rewriting these equations in non-homogeneous variables. For doing this, 
make the change of "time" variable dt' = fj,(x)dt, where fj,(x) is a homogeneous polynomial of degree 
s — 1. Now the system is: 

dxi = A(x) 
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Both the sides of equation (64) are homogeneous expressions on x with the same degree of homogeneity, 
i. e. if x is scaled x' = kx, both the sides of equation will be scaled equall}|^J Denoting £j = 



dt' 



XiXi - xiXi Ai(x)x\ - Ai(x)xi Ai{£) - A 1 (£)£ i 



(xi) 2 



(65) 



The equalities £1 = 1 and £i = were used. In s ]3.7| is proven a statement, that a vector satisfies: 

Ai(x)xj — XiAj(x) = (66) 



iff it is an eigenvector of A. Stationary points of (65) obey the equations: 

Ai(x)x\ — Ai(x)xi = 



(67) 



These equations are the equations (66) with fixed j = 1. 

Proposition 4. All n<yn 2 ^ equations (66) follow from n — 1 equations (66) with some fixed j. 



The system of equations (67) is the system (66) with fixed j = 1. For any other j the proof is 
similar. Let x obey (67). One can consider the solution of the system (64), which at t = coincides 

with x. Then for any i (fj) = 0, see (65). So for any i,j: J, (&) = (f^/g) = 0, and then, 
looking again at (65), Ai(x)xj — XiAj(x) = 0. From (65) it is easily seen that unit map does not move 
points in projective space at all, because for unit map Ai(x)xj — XiAj(x) = and therefore = 0. 
The system (65) can be simplified by choosing ^(x) = (xi) s_1 . Then //(£) = 1 and the system becomes: 



dt' 



(68) 



4.1.1 



n 



2: the equations are reduced to quadratures 



In the case n = 2, there is only one non-homogeneous variable: £ = The system (68 



becomes: 



dt' 



A 2 (0-A 1 (0^ = P(0, 



(69) 



where P(£) is some polynomial of £ of degree s + 1 (or lower, if some coefficients of A vanish). This 
equation is then solved by separation of variables: 



W) 



t' + C 



When £(t') is known, one can find x\(t') by solving: 

dxx A 1 (x) (xi)M!(£) 



dt' 



n(x) 



dxi 



-i 



= *iA 1 (0 
dt' 



(70) 

(71) 
(72) 



x\{t') can be found from this formula. Then: 



t + C 



dt' = {xx{t')) s - l dt 
dt' 



(xi(t')Y 



From these integrals t'{t) can be found, and, recalling that x 2 = £xi, one can return to initial variables 
t,xi(t),X2(t). This case can be used to test some statement and predictions of general theory. 



2 This explains the choose of degree of \x 
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4.2 The "eigenvector" solution 

As it was shown in previous subsection, if Xi(t) is a solution of 

Xi — A^ Xi-^X^Xi^ ' ' ' Xi s (73) 

with initial condition xf , proportional to an eigenvector of A, then Xi(t) evolves so, that 

Xi(t) 



Xl(t) 



const (74) 



This means, that the solution has a form Xi{t) = g{t)xf^ = /(i)e,. If e\ is zero eigenvector, then 
xi = x\ and does not change with time. If ej is unitary eigenvector, then, substituting Xi(t) = f(t)ei 



in (73), one obtains the equation on f(t) with initial condition. In linear and non-linear case the 
equations are different. In linear case: 

f = ww 

f(t)t=o = fo 

Here Aj is eigenvalue, which corresponds to this eigenvector. Solution of this equation: 

/(*) = /oe Alt 



In non-linear case: 



Solution: 



^ = /(*)' (75) 
/(*)t=o = /o (76) 



/(*) = 1 n~r = t^t (77) 



(1 - (* - lX/o)- 1 *)^) 



v/ ( r 

From this formula it is seen, that in non-linear case each solution has a singularity at the moment of 
time to = - — 7vjs=i , i. e. the solution grows infinitely during a finite time. Any stationary point, where 



the linear term vanishes (see 1.4.1), is unstable in its complex vicinity. To show this it is sufficient to 
take as initial condition such a vector, that is proportional to unitary eigenvector and fo > 0. The 
vector of initial condition can be made arbitrary small. Then the solution with this initial condition 
will go at infinity at some finite moment of time. If the map, corresponding to an equation, has a real 
eigenvector, then this point is unstable in mentioned above sense in its real vicinity. 

4.3 Non-linear condition of instability 



Any stable point, in which linear term equals zero (see sect 1.4.1 ) and main non-linear term has unitary 
eigenvector, is unstable in its complex vicinity. To show this, it is sufficient to take as initial condition 
a vector, which is proportional to a unitary eigenvector and fo > 0. This vector can be done arbitrarily 
small. Then the solution will go at infinity in a some finite time. If the map has real unitary eigenvector, 
this point is unstable in its real vicinity. The question of stability of the point of equilibrium of the 
map without unitary eigenvectors is open. This case is encountered rarely. The map should be very 
degenerate, if it has no unitary eigenvectors. For example, the map 

Xl \ / 
x 2 J \ ax\x 2 

has no unitary eigenvectors, and two zero eigenvectors (0, 1) and (1,0). This is the result of solving 



of equations (40) and (39). This map has only one non-zero component of six possible. So, here is a 
brief summary: 
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1. If a map has an unitary eigenvector, the origin is unstable. 

2. If a map has non-zero resultant, the origin is unstable 

3. If a map has no unitary eigenvectors, one cannot make any definite prediction about stabil- 
ity/instability of the origin. 

If the resultant of the map does not equal to 0, the map cannot have zero eigenvectors. So, it has 
unitary eigenvectors and one comes to item 1. 

4.4 The example of equations with vanishing linear term 

Consider the system of differential equations of chemical kinetics of this system of reactions and 
reagents: 

C + D -» A + B 
C + A^ B + D + D 

A, B, C, D represent the reagents. X\, X 2 , X3, X A are, respectively, partial concentrations of -A, B, C, D. 
The condition of mass conservation in this system looks like thi^j 



X 1 + X 2 + X 3 + X 4 = 1 
The equations of chemical kinetics for this system: 

Xi = If 34X3X4 — X31X1X3 + .K24X2 (X4) 2 

X 2 = X34X3X4 + ET31X1X3 — X24X2(X4) 2 
X3 = — X34X3X4 — X31X1X3 + X 2 4X 2 (X 4 ) 2 

X 4 = -X 34 X 3 X 4 + 2X3iX!X 3 - 2^24X2 (X 4 ) 2 , 



(78) 



(79) 



where dot denotes time derivative, #34, #31, #24 are rate coefficients of the reactions. Actually, they 
depend on the temperature and pressure. If these dependencies are neglected, the equations of type 
(63) are obtained. Prom (78) X4 can be expressed through Xi,X2,X 3 , and there will be only three 
variables Xi, X 2 , X3, One of stationary points of this system is: 



X 
X 

x. 



(o) 
1 

(0) 

2 

(0) 



(80) 



Denote now 81 



X 1 -x[ 0) ,6 2 = X 2 -Xi°\5 3 



X3 — Xg \ The linear term in the expansion of 
right-hand sides of (1791) vanishes near the point (1801). In the new variables, 61,62,63, the system (79) 
is: 



S 1 = K 2i 5\ + X 24 <5 2 2 + {K 2i - K u )6l + 2X24M2 + (2X 24 - K M )5 2 5 3 + (2X 24 - X 34 - K 31 )6t6 3 + 

+cubic terms 

5 2 = -K 2A 5l - K 2A 5 2 2 - (X 24 + X 34 )<5 3 2 - 2X24M2 + (#31 - #34 - 2X24)^1 - (2# 24 + #34)^3 + 

+cubic terms 

63 = K 2A 5\ + 2X24^2 + (#34 + 2X 24 - #31)^1 + #24^1 + (2X 24 + K 3A )6 2 6 3 + (# 24 + # 34 )<5| + 

+cubic terms 



;82) 
;ss) 



To analyse the stability of stationary point (80 ), one should consider only the term in expansion of the 
lowest degree, in this case, of second degree. After neglecting cubical terms in (81)-(83), one obtains 
the differential equation with homogeneous quadratic right-hand side, i. e. the equation of type (63). 



3 The mass conservation law is not always so simple. It depends on the particular substances. It can occur that there 
are several independent equations of mass conservation, for each element. Very simplified approximation is considered 
here. 
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4.4.1 Investigation of stability of this example 



To know whether the origin of the system of equations (81 83) is stable or not, one should consider 
the system of equations: 

Si\ ( K 2A b\ + K 2A 5 2 + (K 2A - K u )8l + 2K 2A b x b 2 + (2K 2i - K 3A )6 2 6 3 + (2K 24 - K 3A - K 3l )5 x 5 3 
to - -K u 5\ - K 2A 5 2 - (K 2A + K 3A )5 2 - 2K 24 S 1 S 2 + (K 31 - K 3A - 2X^)5^ - (2K 2A + K 3A )5 2 6 3 | (84) 
S3 J \ K 2A 5\ + 2K 2A 8 X 8 2 + (^34 + 2K 2A - #31)^1 + K 2A 5 2 + (2if 24 + K-^hh + (^24 + K 3A )8 2 

All unitary eigenvectors satisfy the equation: 

<*i = K 2A 5 2 + K 2A 5 2 + (K 2A - K M )5l + 2K 2A 8 X 8 2 + (2K 24 - K u )5 2 5 3 + 

+(2K 24 - K 3A - K 3 i)6 1 5 3 

5 2 = -K 2A b\ - K 2A b\ - (K 2A + K 3A )5 2 3 - 2K 2A b x b 2 + {K 3X - K 3A - 2K 24 )<5 3 <5i - 

(2K 2A + K 3A )5 2 5 3 

8 3 = K 2A b\ + 2K 2A b X b 2 + {K 3A + 2K 24 " K 31 }8 3 8 X + K 2A b\ + (2K 24 + K 3A )b 2 b 3 + 

+ (^24 + K 3A )b\ (85) 

One can notice, that f 2 + f 3 = 0, so the resultant of this map equals 0. So one cannot use the second 



item from sect. 4.3 f 2 + f 3 = 0, and therefore, £3 = —5 2 . Substitute it in first two equations. Then 



one can exclude bi from the equations and obtain the equation with single variable 5 2 : 

2{K 2 A + K 3l K 3A )b 2 2 + {K 2A - 3K 3A - K 3l )b 2 + 1 = (86) 

This equation always has roots, so the point (0, 1,0) is unstable in its complex vicinity, since there are 
complex unitary eigenvectors. It has real roots, when (K 2A — 3K 3A — K 3 i) 2 — 8(iT| 4 + K 3 1K 3 4) > 0, 
and then it is unstable in its real vicinity. 

5 Example of finding degeneracies and peculiarities of non-linear map 



In this section an example of quadratic map of two variables in canonical form is considered. Never- 
theless, all results, formulated in sections |5.1.1[ |5.1.2| |5.1.3[ can be obtained for any non-linear map in 



any coordinates. The only difference (between this example and general situation) is that the devel- 
oped tools are sufficient to fully investigate and classify this example, but in general situation higher 
resultants, higher complanarts, etc can be needed. Only the case when two unitary eigenvectors are 
chosen to be basis vectors is considered. 
Our example is: 

xi -> x\ + 2axxx 2 . . 

x 2 — ► x 2 + 2bx\x 2 

a and b are the parameters, defining a map instead of eigenvalues (see 3.5). 
5.1 General analysis 



By construction of canonical form of the map, (87) has two unitary eigenvectors, namely = (1,0) 
and e^ 2 * 1 = (0, 1). They are unitary eigenvectors under any values of a a nd b. In g eneral case this map 
has c n | s = ^£=p = 3 (s = 2, n = 2) eigenvectors. Here are the equations (39 40|44 ) in this case. 



x\ + 2ax\x 2 = xi . . 

2 . 11™ V^) 



xi + 2bx\x 2 = x 2 
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This is a system (39 ) in this case. One should solve this system of equations to find unitary eigenvectors. 



x\ + 2ax\X2 
x\ + 2bx\X2 








(89) 



This is a system (|40j) in this case. One should solve this system of equations to find zero eigenvectors. 

(90) 



x\ + 2ax\X2 



x\ + 2bx\X2 



X2V, 



This is a system (44) in this case, and all eigenvectors (zero, non-zero unitary) are its solutions. 



Additional homogenizing variable is also called y here, as in (j44J) . A solution of ( |90| ), {x\,X2,y) with 

y 



corresponds to zero eigenvector. This zero eigenvector is (x\,X2)- A solution of ( |90| ), (xi,X2,y) 
with y / corresponds to non-zero eigenvector. This non-zero eigenvector is (xi, X2). Since the system 



(90) is a system of homogeneous equations, the solutions of this system can be scaled by any non-zero 



number. To get unitary eigenvector from a non-zero eigenvector, one should rescale a solution of (90) 



with y / by dividing it by y: then a solution of (90) becomes (xi/y, X2/y, 1) and unitary eigenvector 
is (x 1 /y,x 2 /y). 

5.1.1 Preliminary information about eigenvectors from complanart 



(90) is a system of two homogeneous equations of three variables. Such system possesses a complanart. 



Complanart of this system is calculated in sect |2.3.5[ and it equals 

C = (1 -2a) 4 (l -2b) 4 



When it vanishes, the system (90) has three complanar in the space x\,X2,y roots. A system (90) 



is a system of two quadratic equations, so in general case it has 2-2 = 4 solutions. Two solutions 
of (90) correspond to known unitary eigenvectors, = (1,0) and = (0,1). These solutions 
of (90) are A^ = (1,0,1) and A( 2 ) = (0,1,1). There should be one more eigenvector in this case, 
since the map (87) in general case has c n \ s = c 2 | 2 = 3 eigenvectors. Denote this eigenvector by 

e^ 3 ) = (e^\ e^)- Then the solution of (90 ), corresponding to this eigenvector, is A^ 3 ^ = (e[ , e 2 , e^). 



If is zero eigenvector, ef^ = 0, if e^ 3 ** is non-zero eigenvector, eg 7^ 0, and ef^ = 1 if is 



unitary eigenvector. And the last solution of (90) is A^ 4 ) = (0,0,1), this is a solution (47) for this 



system. Vanishing of complanart means that ^iwsA^A^Ag 3 ) = with at least one triple of indices 
31,32,33, ji 7^ 32, ji / 33,32 / js> Particularly, vanishing of complanart means vanishing of at least 
one of the following expressions: 



Hi 2 i3 A (i) A (2) A (3) _ c (3) _ o (3) 



12 «3 



1 



£ ili2l3 a( 1 ) A ( 2 ) A ( 4 ) 

h 'h h 

r iii 2 i3 A ( 2 ) A (3) A (4) = 

: n^*3 A W A (3) A J4) = 

l\ 12 13 



J 3 ) 



,(3) 
'1 

,(3) 



(91) 
(92) 
(93) 
(94) 



The solution A^ 4 ) enters in (92), (93), (94), so these equations are nothing but conditions of pairs of 



eigenvectors to be complanar, see discussion in sect 3.6 The equation (92) controls complanarity of 



and e@\ (93) - of and and (94) - of eW and e^. The condition (91) has no such a simple 
interpretation. So, if a = 1/2 or b = 1/2, complanart vanishes, and one expects that either and 
or and e^ 3 ) will be collinear, or the condition (91) will hold. 



5.1.2 Preliminary information about eigenvectors from resultant 

R{A} = 1 - Aab 



Resultant of the map (87) equals: 
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As it was stated in sect |3.4| zero eigenvectors exist iff the resultant of the map equals 0. So one expects, 
that if 1 — 4a6 = 0, there will be zero eigenvectors, and if 1 — 4a6 / 0, there will be no zero eigenvectors, 
and will be only non-zero eigenvectors. 



5.1.3 Preliminary information about eigenvectors from the equation XiAj —xjA{ 



As it was stated in sect 3.7.1 if polynomials XiAj(x) — XjAi(x) are identically zero, the map A is unit 



map. For the map (87), there is only one equation: 

x\A2(x) — X2A\(x) = x\{x\ + 26x1x2) — xi(x\ + 2ax\X2) = x\x\(\ — 2a) — x\y(l — 2b) = (95) 

So, if a = 1/2 and 6 = 1/2, the map A will be unit map, every vector of the space will be eigenvector 
of A, every polynomial of degree s — 1 fi(x) will be an eigenvalue of the map A, and characteristic 
polynomial of A will be identically zero. 



5.2 Eigenvectors 



Besides two solutions (xi,x%,y) = (0,1,1) and (xi,X2) = (1,0,1), the system (90) has the following 
solution: 



(96) 



X\ \ 




( \ 


- 2a 




H 




- 26 


y J 






-4a6 



If y value is equal to zero, (96) corresponds to zero eigenvector: 

V W = ( !- 2 « 
~ V 1 - 26 

(3) 

(of course, the normalization is arbitrary). In this case = 0. 



(97) 



If y value is not equal to zero, (96) corresponds to non-zero eigenvector. In this case, the solution 
vector can be renormalized: 



l-2a 
l-4a& 
l-2b 
1-Aab 
1 



(98) 



to get unitary eigenvector: 



,(3) 



l-2a 
1-Aab 
1-26 
l-Aab 



In this case = 1. In agreement with predictions of sect 5. 1.2[ zero eigenvector exists iff 1 
or when the resultant vanishes. In agreement with predictions of sect 5.1.1 if a = 1/2 or 6 
are two coinciding eigenvectors. If a = 1/2, (99) coincides with 
with 



The prediction of sect 



5.1.3 



and if 6 = 1/2, (99) coincides 
namely, the case a = 1/2 and 6 = 1/2 requires separate 



(99) 

-4a6 = 0, 
1/2, there 



consideration. All the components of (96) are now zero. The equation on eigenvectors/eigenvalues 
with a = 1/2,6 = 1/2: 



Xf + X]X2 
2 + XiX 2 



X 



xiy 
x 2 y 



(100) 



The case of unit map is the case of decreasing range of the system (90). Unfortunately, there have 
been developed no clear methods for determining whether the system of non-linear equations is degen- 
erate. Complanart is an attempt to full this gap, but it is easily seen from this example that it does 
not distinguish the cases when two vectors coincide and the decreasing of range of a system. Only 
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eigenvectors, which are not equal to known (0, 1) and (1,0), are interested in now. So, first equation 
of (100) can be divided by x\ and the second one can be divided by x 2 . 



xi + x 2 = y 
xi + x 2 = y 

The equation x\ + x 2 = y has projective one-dimensional space of solutions: one can set y = 1 (to 
obtain unitary eigenvectors), x 2 = C and x\ will be 1 — C. For any C the vector 

1 -c ) (101 » 

is unitary eigenvector. But there is also zero eigenvector, (xi,x 2 ,y) = (—1,1,0). Arbitrary vector of 
the space (x\,x 2 ) is zero eigenvector, if x\+x 2 = 0, and (xi, x 2 ) is non-zero eigenvector, if x\ + x 2 ^ 0. 

5.3 Eigenvalues 

Pol ynomial A, corresponding to unitary eigenvector eW is found from equation A(e^) = 1, see sect 
3.2 This equation for e^ 1 ) is Ai = 1, i. e. to the eigenvector e^ 1 ) corresponds the line Ai = 1 on the 
plane Ai,A2- And for e^ 2 ': A2 = 1, i. e. to the eigenvector corresponds the line A2 = 1 on the 
plane Ai, A 2 . For e^: X 1 ^ B + ^A 2 = 1, or 

Ai(l - 2a) + A 2 (l - 26) = 1 - 4a6 (102) 

This is also a line at Ai,A2 plane. To find A, corresponding to zero eigenvector, one should solve 



A(t;W) = 0, i. e. in this case Ai(l — 2a) + A2(l — 26) = 0. But (102) is equivalent to this equation, 
because the resultant, namely 1 — 4a6, turns to zero if there is a zero eigenvector. It is not a random 
coincidence, namely the term Ai(l — 2a) + A2(l — 26) — (1 — 4a6) is a factor in the decomposition of 
characteristic polynomial. It is easily seen, that in the case of a = 1/2 and 6 = 1/2 (or when our map 
is unit map), any A is eigenvalue. 

5.4 Characteristic polynomial 



The system of equations (37) for the map (87): 



x\ + 2ax\x 2 = (A1X1 + \ 2 x 2 )x\ 
x\ + 2bx\x 2 = [\\x\ + A 2 x 2 )x 2 

Characteristic polynomial is the resultant of the system: 

x\ + 2axiX2 — (A1X1 + A2X2)xi = 

x\ + 26xiX2 — (A1X1 + A2X2)X2 = 

The characteristic polynomial is equal: 

Ch A {\) = (1 - Ai)(l - A 2 )(l - 4a6 - Ai(l - 2a) - A 2 (l - 26)) (103) 

Here we would like to emphasize one more time: decomposability of characteristic polynomial is non- 
trivial property. The polynomial 1 + (Ai) 2 + (A2) 2 , for example, cannot be decomposed on linear on 
Ai, A2 factors. If the resultant does not equal to 0: 

Ch A (X) = (1 - 4a6)(l - A(e (1) ))(l - A(e (2) ))(l - A(e< 3 >)) 



in full agreement with (46). If resultant of a map is equal to zero: 

Ch A {\) = (1 - A(e«))(l - X(e^))X(v^) 



in full agreement with (43). If a = 1/2 and 6 = 1/2 (i. e. unit map), characteristic polynomial is 



identically zero, in full agreement with s 3.7.1 
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5.5 Phase diagram 



Differential equation, which corresponds to the map (87) , is: 



x\ = x\ + 2ax\X2 

x\ = x\ + 1bx\X2 



(104) 



5.5.1 a / 1/2, b / 1/2,1 - 4ab ^ 

e of any 
)• Figgi 



This is the case of absence of any degeneracies and pecularities. There are three unitary eigenvectors, 
(1,0); (0, 1); ( \zlab ' i-Iab ) 1 Figjl| is the phase portrait of this system with a = —1, b 



■1. 




Figure 1: The phase portrait with a = b = — 1. There are no degeneracies and peculiarities, there are 
three unitary eigenvectors: (0, 1), (1, 0), (3/5, 3/5) 



5.5.2 a/ 1/2, 1/2,1 - 4ab = 

The map in this case is degenerate, i. e. it has a zero eigenvector (1 



with sect 



5.1.2 



Each point, which lie on the line j 



-2a 



2a, 1 — 26) with full accordance 
Yzjzb {^is is a line of zero eigenvectors), is a 



._ 



stationary point. Figj2]is the phase portrait of this system with a = —1, b = —1/4. 
5.5.3 a/l/2,b = l/2 

Because the complanart of this map equals 0, there is a double eigenvector of this map, with full 



accordance with sect 5.1.1 In this case the map has only two unitary eigenvectors, (0, 1) is "simple 1 



eigenvector (with multiplicity 1), and (1,0) with multiplicity 2. The eigenvectors with multiplicity 2 
have a special property: the phase trajectories tend in projective space to this eigenvector from one 
side, and they tend out of it in projective space at other side of this eigenvector. On figj3]at upper right 
corner trajectories tend to eigenvector (1,0) in projective space (namely, xi — ► 00,22 — * 00, — > 0) 
and at lower right corner trajectories move out of eigenvector (1,0) in projective space. The equation 
on X2/X1 = C in the vicinity of looks like 



(105) 
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Figure 2: The phase portrait with a = — 1, b = —1/4. The map is degenerate, so it has zero eigenvector 
(3,2). All the points of the line § = § are stationary points. 



where C is a constant, depending on a and b. This projective equation explains such a behaviour: ( 
has the same signs at both sides from stationary point, so the solution approaches from one side and 
move away from other side. Near a double eigenvectors, the linear term in projective equation always 
vanishes. Figj3]is the phase portrait of the system with a = 1, b = 1/2. The case with a = 1/2, b ^ 1/2 
reduces to this case by substitution x\ — * X2, X2 — ► x\. 



5.5.4 a = l/2,b = l/2 

According to predictions of 5.1.3 this map is unit map and any vector is an eigenvector. Figj4] is the 
phase portrait in this case. 



5.5.5 Phase diagram without unitary eigenvectors 

The example of map without unitary eigenvectors is considered here. There are no methods to deter- 
mine whether the point which has not unitary eigenvectors is stable or not. The example is: 

Xl ^ , ° (106) 

X2 — > bX\X2 

The basis consists of two zero eigenvectors. Besides these vectors, there are no other eigenvectors. In 
this particular case the point (0, 0) is unstable. Phase diagram is at the Figj5} 
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Figure 3: The phase portrait with a = 1, b = 1/2. Because the complanart of the map equals 0, there 
is double eigenvector (1,0). Trajectories at upper right corner tend to eigenvector (1,0) in projective 
space and at lower right corner trajectories move out of eigenvector (1,0) in projective space. This is 
a general property of double eigenvectors. 




Figure 4: The phase portrait with a = 1/2,6 = 1/2. Unlike all previous cases, all the phase tra- 
jectories are straight lines. This is because this map is unit map, and all vectors of the space are 
eigenvectors. Vectors, proportional to (—1,1) vector, are zero eigenvectors, and all other vectors are 
non-zero eigenvectors. 
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